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Abstract 


We calcnlate qnark and lepton masses and qnark mixing angles in the framework of 
a snpersymmetric SU(4)®SU(2)^®SU(2)j:j model where the gange gronp is broken 
at 10^® GeV. The model predicts third family top-bottom-tan Ynkawa nnihcation, as 
in minimal SO(IO). The other smaller Yukawa couplings are assumed to arise from 
non-renormalisable operators suppressed by powers of some heavy scale. We perform 
a systematic operator analysis of the model in order to hnd the minimum set of 
operators which describe the low energy quark and lepton masses, and quark mixing 
angles consistent with low-energy phenomenology. A novel feature of the model is 
the possibility of asymmetric texture zeroes in the Yukawa matrices at the scale of 
the new physics. Successful predictions are obtained for m*, tan/3, m^/me 

and Vub in terms of a CP violating phase cf. For example, we predict tan j3 = ?>b — 65, 
rut = 130 — 190 GeV, and |Wb| > 0.0040. 


1 Introduction 


The problem of quark and lepton masses is one of the most fascinating and perplexing 
problems of particle physics. The standard model, despite its successes, can offer no 
glimpse of insight into the apparently bizarre pattern of masses and mixing angles 
which experiment has presented us with. We do not even know why there are three 
families rather than one. It is clear, then, that in order to gain some insight into the 
fermion mass problem, one must go beyond the standard model. The big question of 
course is what lies beyond it? 

We have not yet experimentally studied the mechanism of electroweak symmetry 
breaking, so one might argue that it is premature to study the fermion mass problem. 
Unless we can answer this, we have no hope of understanding anything about fermion 
masses since we do not have a starting point from which to analyse the problem. 
However LEP has taught us that whatever breaks electroweak symmetry must do 
so in a way which very closely resembles the standard model. This observation by 
itself is enough to disfavour many dynamical models involving large numbers of new 
fermions. By contrast the minimal supersymmetric standard model (MSSM) mimics 
the standard model very closely. Furthermore, by accurately measuring the strong 
coupling constant, LEP has shown that the gauge couplings of the MSSM merge 
very accurately at a scale just above 10^® GeV, thus providing a hint for possible 
unihcation at this scale. On the theoretical side, supersymmetry (SUSY) and grand 
unified theories (GUTs) fit together very nicely in several ways, providing a solution 
to the technical hierarchy problem for example. When SUSY GUTs are extended 
to supergravity (SUGRA) the beautiful picture of universal soft SUSY breaking pa¬ 
rameters and radiative electroweak symmetry breaking via a large top quark yukawa 
coupling emerges. Finally, there is an on-going effort to embed all of this structure 
in superstring models, thereby allowing a complete unihcation, including gravity. 

Given the promising scenario mentioned above, it is hardly surprising that many 
authors have turned to the SUSY GUT framework as a springboard from which to 
attack the problem of fermion masses [Q]. Indeed in recent years there has been a 
hood of papers on fermion masses in SUSY GUTs. Although the approaches diher 
in detail, there are some common successful themes which have been known for some 
time. For example the idea of bottom-tan Yukawa unihcation in SUSY GUTs P| 
works well with current data [Q. A more ambitious extension of this idea is the 
Georgi-Jarlskog (GJ) ansatz which provides a successful description of all down-type 
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161, and which also works well with current 


quark and charged lepton masses 
data 0. The GJ approach involves the idea of texture zeroes and predicts simple 
relations at the unihcation scale which are then evolved to low-energies using the 
renormalisation group (RG) equations. These approaches are concerned with general 
properties of the mass matrices, rather than those of specihc models. 


In order to understand the origin of the texture zeroes, one must consider the 
details of the model above the scale Mx 10^® GeV (SO (10) for example in the 
case of GJ). While one might not wish to restrict oneself to some particular gauge 
symmetry at Mx, it is almost essential to specify the model at this scale in order to 
make any progress at all. The alternative is to simply make a list of assumptions about 
the nature of the Yukawa matrices at Mx 0. For example Ramond, Roberts and 
Ross (RRR) [0 assumed symmetric Yukawa matrices at Mx, together with the GJ 
ansatz for the lepton sector. It is difficult to proceed beyond this without specifying 
a particular model. Indeed, this model dependence may be a good thing since it may 
mean that the fermion mass spectrum at low energies is sensitive to the theory at 
Mx, so it can be used as a window into the high-energy theory. Therefore in what 
follows we shall restrict ourselves to the very specihc gauge group at Mx referred to 
in the title. Our motivation for considering this particular theory is discussed below. 


Twenty-one years ago Pati and Salam proposed a model in which the standard 
model was embedded in the gauge group SU(4)®SU(2)2,®SU(2)ij 0. More recently 
a superpersymmetric (SUSY) version of this model was proposed in which the gauge 
group is broken at Mx 10^® GeV 0. The model 0 does not involve adjoint 
representations and later some attempt was made to derive it from four-dimensional 
strings, although there are some difficulties with the current formulation |^. In this 
paper we shall not be concerned with the superstring formulation of the model, but 
instead we shall focus on the “low-energy” effective held theory. The absence of ad¬ 
joint representations is not an essential prerequisite for the model to descend from the 
superstring, but it leads to some technical simplihcations. Also in the present model, 
the colour triplets which are in separate representations from the Higgs doublets, be¬ 
come heavy in a very simple way so the Higgs doublet-triplet splitting problem does 
not arise. These two features (absence of adjoint representations and absence of the 
doublet-triplet splitting problem,) are shared by hipped SU(5)(8)U(1) 0, which also 
has a superstring formulation. Although the present model and hipped SU(5)G)U(1) 
are similar in many ways, there are some important diherences. Whereas the Yukawa 
matrices of hipped SU(5)®U(1) are completely unrelated at the level of the ehective 
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field theory at Mx (although they may have relations coming from the string model) 
in the present model there is a constraint that the top, bottom and tan Yukawa 
couplings must all unify at that scale. In addition there will be Clebsch relations 
between the other elements of the Yukawa matrices, assuming they are described by 
non-renormalisable operators, which would not be present in flipped SU(5)®U(1). In 
these respects the model resembles the SO(IO) model recently analysed by Anderson 
et al |T^. However it differs from the SO(IO) model in that the present model does 
not have an SU(5) subgroup which is central to the analysis of the SO(IO) model. 
In addition the operator structure of the present model is totally different. Thus 
the model under consideration is in some sense similar to flipped SU(5)®U(1), but 
has third family Yukawa unihcation and precise Clebsch relationships as in SO(IO). 
We hnd this combination of features quite remarkable, and it seems to us that this 
provides a rather strong motivation to study the problem of fermion masses in this 
model. 

The problem of fermion masses in the SU(4)®SU(2)j;,(8)SU(2)j:j model was recently 
considered by one of us |TB]. It was implicitly assumed that the Yukawa matrices were 
symmetric, and it was shown that by introducing suitable operators the model could 
make contact with the successful RRR ansatz, which incorporate the GJ ansatz, at 
the expense of hne-tuning the coefficients of the operators El- Essentially the small 
entries in the RRR matrices were obtained by assuming that the coefficients of two op¬ 
erators were tuned to partially cancel. The purpose of the present paper is threefold. 
First we shall generalise the analysis to the case of non-symmetric Yukawa matrices, 
since there is no symmetry which enforces symmetric Yukawa matrices in this modelQ. 
This allows the possibility of asymmetric texture zeroes, which as far as we are aware 
have never been considered before. Of course this means that we cannot rely on the 
RRR analysis, and therefore we perform our own phenomenological analysis of the 
quark and lepton masses and quark mixing angles. Second we extend the operator 
analysis to consider many other operators not considered in the previous analysis. In 
fact we search for all possible low dimensional operators, and systematically search 
for the minimum set with which to describe the spectrum. Third we impose a nat¬ 
uralness criterion and reject all possibilities which involve fine-tuning the coefficients 
of operators. The result of all this is a small set of possible solutions to the problem 
of quark and lepton masses and CKM angles in this model. 

The remainder of the paper is organised as follows. In section 2 we briefly sum- 


'^Even the imposition of parity does not lead to symmetric Yukawa matrices (see Section 3.3.) 
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marise the model. In section 3 we describe the operator strategy we employ. In 
section 4 the details of the calculation are outlined, including the RG and CKM 
analysis. The ansatze, results and predictions are presented. Section 5 contains our 
conclusions about the previous analysis, and a brief discussion of theoretical uncer¬ 
tainties involved with the calculation. In the appendices we list the operators in 
explicit component form. 


2 The Model 


Here we briefly summarise the parts of the model which are relevant for our analysis. 
For a more complete discussion see 0. The gauge group is, 

SU(4) ® SU(2)^ ® SU(2)^. (1) 

The left-handed quarks and leptons are accommodated in the following representa¬ 
tions, 

/ p p \ i 

_ / > 11 ^ 1 !^ 1 !^ 7 / \ 

( 2 ) 
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(3) 

where a = 1... 4 is an SU(4) index, a,x = 1, 2 are SU(2)x, indices, and i = 1... 3 
is a family index. The Higgs helds are contained in the following representations. 




(4) 


(where hi and ^2 are the low energy Higgs superhelds associated with the MSSM.) 
The two heavy Higgs representations are 


A, = (4,1,2)= (ft ft f| f« 


ug uf ug Uh 


iH 


and 


j 

The Higgs helds are assumed to develop VEVs, 

< H >=< Ufj >~ Mx, < H >=< Pjj >~ Mx 
leading to the symmetry breaking at Mx 

SU(4) ® SU(2)^ ® SU(2)^ —. SU(3)^ (g) SU(2)^ ® U(l)y 


(5) 


( 6 ) 


(7) 


( 8 ) 
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in the usual notation. Under the symmetry breaking in Eq.§, the Higgs held h in 
Eq.^ splits into two Higgs doublets hi, h 2 whose neutral components subsequently 
develop weak scale VEVs, 


< h\ >= Ui, < ^2 >= V2 


(9) 


with tan/? = V 2 /V 1 . 

In addition to the Higgs helds in Eqs. the model also involves an SU(4) sextet 
held D = (6,1,1) and three singlet helds 0 ^ = ( 1 , 1 , 1 ) which do not acquire VEVs 
plus one singlet held N = (1,1,1) which acquires a weak scale VEV < N >= x. The 
superpotential, suppressing gauge indices, is then 

W = Ai^FiT;h+ArEiih0^+A3MT) + A4MT> + A5hhiV+A;r''Vm0n05+A7iV3 (10) 


Note that this is not the most general superpotential invariant under the gauge sym¬ 
metry. Additional terms not included in Eq.TD may be forbidden by imposing suitable 
discrete symmetries, the details of which need not concern us here. The D held does 
not develop a VEV but the terms in Eq. 0 HHD and HHD combine the colour 
triplet parts of H, H and D into acceptable GUT-scale mass terms . The 0^ helds 
play an important part in ensuring that the right-handed neutrinos gain large masses, 
leading to acceptably small observable neutrino masses. The ehect depends on terms 
in the superpotential like Fiif0 and 0^ |T^. Below Mx the part of the superpotential 
involving quark and charged lepton helds is just 


W = >!SQiUjh2 + X%QiDjhi + X%LiEjhi + ... 


( 11 ) 


with the boundary conditions at Mx, 

= 4 = X^ = yi, A 5 = A, Xr = k (12) 


The model just described must explain why the gauge couplings which are roughly 
equal at Mx ~ 10^® GeV remain roughly equal up to the compactihcation scale 
Me ~ 10^^ GeV. Gonventional SUSY GUTs achieve this in the most direct way 
possible, by embedding the standard gauge group into a simple gauge group with 
a single gauge coupling constant. However, conventional SUSY GUTs are not fully 
unihed because they do not include gravity. The only consistent known theories of 

^The resulting low energy theory may resemble the so-called Next-to-Minimal Supersymmetric 
Standard Model (NMSSM) involving an extra gauge singlet. However for simplicity our calculations 
will be based on the MSSM. 
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gravity are string theories, and string theories which allow adjoint superhelds are 
quite cumbersome ||^. On the other hand string theories that do not involve adjoint 
superhelds, and consequently cannot involve a simple gauge group, must explain why 
the gauge couplings which appear to be unihed at Mx are in fact unihed at Me- 


Recently it was suggested by one of us that an attractive solution to this 
problem is to introduce some additional GUT-scale superhelds in order to make the 
model left-right symmetric. 


//'= (4,2,1), //'=(4,2,1) 


(13) 


Having guaranteed the equality of the SU(2)£,^/j couplings ql = Qr, it is possible to 
require that the SU(4) beta function, /34, is equal to the common SU(2) beta functions, 
P 2 , so that if the gauge couplings are equal at Mx then they will remain equal above 
this scale. The one-loop SUSY (3 functions are. 


A 


dai 
^ dii 



(14) 


where in the model dehned in the previous section, and augmented by the Higgs 
representations in Eq.|T^ we hnd 


bt = [6 - no - ino], 62 = [-1 - inn] 


(15) 


where we have allowed for ud copies of the sextet superheld H, and nn copies of 
the set of helds [H, H, H', H'). From Eqs. |T^, ^ it is clear that the combination 
of left-right symmetry and the choice ud = 7 (for any choice of rin) is sufficient to 
guarantee that if the couplings are equal at Mx then they will remain equal above 
this scale to one-loop order, ignoring threshold effects. However, as we shall see in 
section such a left-right symmetry does not lead to any simpliheations of the 
Yukawa matrices at Mx, and so we shall not impose such a symmetry in this paper. 


3 Operators 

3.1 The Basic Strategy 

In this model the two Higgs doublets are unihed into a single representation h in Eq.^ 
and this leads to the GUT-scale equality of the three Yukawa matrices in Eqs.pdj, pYj 
This boundary condition also applies to the version of the conventional SUSY GUT 
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based on SO(IO) in which both Higgs doublets are unihed into a single 10 represen¬ 
tation. As it turns out, the idea of Yukawa unihcation works rather well for the third 
family |^, leading to the prediction of a large top quark mass rrit > 165 GeV, and 
tan/3 ~ mt/rrih where mf, is the bottom quark mass. However Yukawa unihcation 
for the hrst two families is not successful, since it would lead to unacceptable mass 
relations amongst the lighter fermions, and zero mixing angles at Mx- In the SO(IO) 
SUSY GUT there are various ways out of these difficulties, and if the present model 
is to be regarded as a surrogate SUSY GUT it must also resolve them. 


One interesting proposal has recently been put forward to account for the fermion 
masses in an SO(IO) SUSY GUT with a single Higgs in the 10 representation 
According this approach, only the third family is allowed to receive mass from the 
renormalisable operators in the superpotential. The remaining masses and mixings 
are generated from a minimal set of just three specially chosen non-renormalisable 
operators whose coefficients are suppressed by some large scale. Furthermore these 
operators are only allowed to contain adjoint 45 Higgs representations, chosen from 
a set of helds denoted 45y, 45b-l, 45x whose VEVs point in the direction of 

the generators specihed by the subscripts, in the notation of [Q. 

This is precisely the strategy we wish to follow here. We shall assume that only the 
third family receives its mass from a renormalisable Yukawa coupling. All the other 
renormalisable Yukawa couplings are set to zero. Then non-renormalisable operators 
are written down which will play the role of small effective Yukawa couplings. The 
effective Yukawa couplings are small because they originate from non-renormalisable 
operators which are suppressed by powers of the heavy scale M. In this paper we 
shall restrict ourselves to all possible non-renormalisable operators which can be con¬ 
structed from different group theoretical contractions of the helds: 
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Oy ~ (F,Fj)h 



+ h.c. 


( 16 ) 


where we have used the helds H,H in Eqs.^J^ and M is the large scale M > Mx- 
The idea is that when H, H develop their VEVs such operators will become ehective 
Yukawa couplings of the form hFF with a small coefficient of order We 

shall only consider up to n = 2 operators here, since as we shall see even at this level 
there are a wealth of possible operators that are encountered. Although we assume no 
intermediate symmetry breaking scale (i.e. SU(4)®SU(2)/,®SU(2)/{ is broken directly 
to the standard model at the scale Mx) we shall allow the possibility that there are 
diherent higher scales M which are relevant in determining the operators. For example 
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one particular contraction of the indices of the fields may be associated with one scale 
M, and a different contraction may be associated with a different scale M' . We shall 
either appeal to this kind of idea in order to account for the various hierarchies present 
in the Yukawa matrices, or to higher dimensional operators which are suppressed by 
a further factor of M. 


3.2 A Simple n = 1 Example 


In the present model, although there are no adjoint representations, there will in gen¬ 
eral be non-renormalisable operators which closely resemble those in SO(10) involving 
adjoint helds. The simplest such operators have already been considered in ref. 


and correspond to n = 1 in Eq. with the {HH) group indices contracted together. 


These operators are similar to those of 
adjoint Higgs representations. It is useful to dehne the following combinations of 
helds, corresponding to the different HH transformation properties under the gauge 
group in Eq.|I|, 


[1^ but with HH playing the role of the 


{HH)a = 

(1,1,1) 


{HH)b = 

(1,1,3) 


{HH)c = 

(15,1,1) 

(17) 

{HH)d = 

(15,1,3) 



The explicit form of the operators is given in Appendix 1. For the adjoint combina¬ 
tions we may write B = B^T^, C = C^T^, D = where {x = 1,..., 3) 

are the generators of SU(2)/^, and {p = 1,, 15) are the generators of SU(4). It 
is clear that when the helds H, H develop the VEVs in Eq.|^ the composite helds in 
Eq.[T^ acquire VEVs 


where 


<HH>a = 

< Vh >< > 

<HH>b = 

< vh >< vh > 

<HH>c = 

< Vh >< 

<HH>d = 

<vh><vh> 

rj^B — L _ 

2v/gdiag(l, 1,1, 3) 

rj^ZR _ 

^diag(l,-l). 


(18) 


( 19 ) 


8 





Armed with the above results it is straightforward to construct the operators of the 
form of Eq.|^ explicitly, and hence deduce the effect of each operator. For example 
for n = 1 the four operators are, respectively. 


O 


A,B,C,D 

ij 


FiFjh 


(HH) 


A,B,C,D 


M2 


+ H.C. 


( 20 ) 


where we have suppressed gauge group indices. When the combinations A, B,C,D 
in Eq.|^ acquire the VEVs in Eq.|^ the generators and in Eq.|^ then 

just count the quantum numbers of the components of the helds F, F, leading to 
quark-lepton and isospin splittings, as shown explicitly below: 


Of'j = aij{QiUjh 2 QiDjhi LiEjhi H.c.) 

Ofj = bij{QiUjh 2 — QiDjhi — LiEjhi -\- H.c.) 

Ofj = Cij{QiUjh 2 -\- QiDjhi — 3LiEjhi -|- H.c.) 

= dij{QiUjh 2 — QiDjhi SLiEjhi H.c.) (21) 

where the coefficients of the operators a^, bij, Cij, dij are all of order • 


3.3 Parity 


In ref. |j^ combinations of the operators in Eq.^ were used to reproduce the successful 
RRR and GJ textures. However it is clear that there is no real justihcation for 
assuming that the Yukawa matrices are symmetric in this model. To illustrate the 
point, let us impose a left-right (parity) symmetry on the model of the kind introduced 
earlier in order to ensure that ql = Qr above Mx- Under the Z 2 parity we have. 


Al - 
F* ^ H 
H ^ H' 

H ^ H' 

h ^ h (22) 

where the helds are the gauge helds of SU(2)2 ,,_r and H',H' are (4,2,1) and 
(4, 2,1) irrep.s respectively. It is clear that operators such as those in Eq.|^ will not 
then lead to symmetric Yukawa matrices since under parity as in Eq.^, the {HH) 
combination which develops the VEV and leads to the ehective Yukawa coupling is 
transformed to {H'H') which cannot attain a VEV if electroweak symmetry is to 
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Operators 

Combination 

0 ^ to 0 ^ 

{FF) (HH) h 

(15 0 1,2,2) (15 0 1,1,10 3) (1,2,2) 

0 ^ to 0 ^ 

(FH) (FH) h 

(6 0 10,2,2) (6 0 10,1,10 3) (1,2,2) 

0 ^ to 0 ^ 

(FH) (FH) h 

(1501,2,2) (1501,1,103) (1,2,2) 

0 ^ to 0 ^ 

Mixed group structure 


Table 1: Operator classification of the n 
(see Appendix 1 for more details). 


1 operators, including those in Eq.|T^,[T7 


remain intact at Mx and so does not lead to an effective Yukawa coupling. Note that 
this argument only applies to the non-renormalisable operators. The renormalisable 
operators would lead to symmetric Yukawa matrices if parity was imposed, since h 
is transformed into itself. It is possible that there may be some non-renormalisable 
operators which would lead to symmetric Yukawa matrices at Mx but these are not 
the kind of operators we consider here. It is clear that the nature of this model is to 


lead to non-symmetric Yukawa matrices, so the analysis of ref. 113] must be extended. 


3.4 General Analysis of n = 1 and n = 2 Operators 

The n = 1 operators are by definition all of those operators which can be constructed 
from the five fields FFhHH by contracting the group indices in all possible ways, 
as discussed in Appendix 1. Here we only summarise the results of this analysis by 
listing the group theoretical contractions of fields in Table |^, and the precise group 
structure of the operators in Table |^. After the Higgs fields H and H develop VEVs at 
Mx of the form = vhi {Sax) = {S 41 ) = the operators listed in the 

appendix yield effective low energy Yukawa couplings with small coefficients of order 
Mx/M"^. However, as in the simple example discussed previously, there will be precise 
Clebsch relations between the coefficients of the various quark and lepton component 
fields. These Clebsch relations are summarised in Table Having discussed the origin 
of the effective Yukawa terms in some detail for the {HS) contracted operators, we 
shall now be more schematic in our description of the remaining types of operator. 
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Operator 

SU(4), 

SU(2)i 

SU(2)^ 

0^ 

1® 1 

2 ® 2 

2 ® 1 ® 2 

0^ 

1® 1 

2 ® 2 

2 ® 3 ® 2 

0^ 

15® 15 

2 ® 2 

2 ® 1 ® 2 

0^ 

15® 15 

2 ® 2 

2 ® 3 ® 2 

qE 

6 ® 6 

2 ® 2 

2 ® 2 ® 1 

0^ 

6 ® 6 

2 ® 2 

2 ® 2 ® 3 

0^ 

10 ® 10 

2 ® 2 

2 ® 2 

0^ 

10 ® 10 

2 ® 2 

2 ® 2 ® 3 

0^ 

1® 1 

2 ® 2 

2 ® 2 

0^ 

1® 1 

2 ® 2 

2 ® 2 ® 3 

qk 

15® 15 

2 ® 2 

2 ® 2 

0^ 

15® 15 

2 ® 2 

2® 2® 3. 


Table 2: The combinations shown in Table ^lead to the group structure for the n = 1 
operators shown here, where the group singlet contraction is taken. 



QUh2 

QDhi 

LEhi 

0^ 

1 

1 

1 

0^ 

1 

-1 

-1 

0^ 

1 

1 

-3 

0^ 

1 

-1 

3 

qE 

0 

1 

0 

0^ 

1 

-4 

0 

0^ 

0 

1 

2 

0^ 

2 

1 

2 

0^ 

0 

0 

0 

0'^ 

0 

0 

1 

qK,N,0 

1 

0 

0 

0^ 

5 

1 

3/4 

qM 

0 

1 

1 

0^ 

4 

4 

3 


Table 3: When the Higgs helds develop their VEVs at Mx, the n = 1 operators lead 
to the effective Yukawa couplings with Clebsch coefficients as shown. These results 
are a generalisation of Eq.^. 
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QUH2 

QDHi 

LEHi 

-Q^ 

1 

3 

9/4 

QDd 

1 

3 

3 

QMd 

1 

3 

6 

0 ^ 

0 

1 

1 

0 ^ 

0 

1 

3/4 

0 ^ 

0 

1 

2 


Table 4: n = 2 operators utilised, where is any one of QJ^^^Dp,Dq,Dr,Ds.^ q 2 jg 
one of and is one of operators are explicitly 

defined in Appendix 2. 

In Table ^ we have neglected terms involving the right handed neutrinos since 
they receive a large mass through the see-saw mechanism]^ Note that associated with 
each operator is only one coupling constant, so that for example gives the same 
Yukawa coupling to the up and down quarks and the charged leptons at M^, as in 
Eq.0. No non-renormalisable terms with HH or HH in the Higgs structure are 
both gauge invariant and give a non zero mass term, so only {HH) operators are 
considered. 

The n = 2 operators are by dehnition all those operators which can be constructed 
from FFhHHHH by contracting the group indices in all possible ways, as discussed 
in Appendix 2. There are 400 n = 2 operators, formed by different combinations of 
the SU(4)c and SU(2)ij structures e.g. we label an operator with structures A and 
t as Brevity prevents us from listing all the possible operators, but the useful 
operators are listed in Table Two features of the n = 2 Clebsch coefficients listed 
are useful: the operators down quark Clebsch coefficient to be three 

times that of the up quark, which helps to predict a small |Vufe|, and the operators 
labeled give masses to down quarks and leptons but not to up quarks, helping 

to account for the up-down mass splitting in the hrst family (see section 4.5). 

4 The Calculation 


^Which explains why ~ 0, it only contributes to a uh term. In this paper we shall not consider 
the problem of neutrino masses. 
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3.4 



Figure 1: Running of Masses Between rrif and rrit (which is displayed in Fig.|[) 


4.1 Masses and Mixing Angles at Low Energies 


To constrain the Yukawa matrices at Mx, we need to use renormalisation group 
equations to evolve low energy parameters such as CKM matrix elements and fermion 
masses up to Mx- We denote running fermion masses at the MS scale /r as mj(/r). 


The masses of the fermions are first run up to the top mass rfit using effective 3 
loop QCD ® 1 loop QED in the MS renormalisation scheme ||T7| , |^. Fig.|^ shows 

the parameter rji dehned by 


Vf 


m/(max(m/, IGeV)) 
rfifirfit) 


(23) 


where max(mj, 1 GeV) is the greater of ruf and 1 GeV. 

For given values of tan/3, as{Mz) and MS fermion masses dehned by ffif = 
m/(max(m/, IGeV)), the diagonal Yukawa couplings at frit are determined by 




(24) 

vr]u,c,t sin (3 


V^rfid,s,b 
vr]d,s,b cos [3 

(25) 



(26) 

vrje,ti,T cos /3 ’ 


where v = vl + v 2 = 246 GeV. All values used for the masses are running values, as 


13 

















Lower Bound/GeV 

Upper Bound/GeV 

rud (1 GeV) 

0.0055 

0.0115 

rhs (1 GeV) 

0.105 

0.230 

ifiu (1 GeV) 

0.0031 

0.0064 

rhc{fhc) 

1.22 

1.32 

rhb{fhb) 

4.1 

4.4 


Table 5: Running masses of the lightest hve quarks as provided by m- 



alpha_s (M_Z) 


Figure 2: Physical mass on the top quark predicted by the renormalisable 33 Yukawa 
term in Eq.^ for various values of ruf,. The pole mass is given by mt(pole) = rht[l + 
Aasifht )/to one loop. 


in ref. The values of CKM mixing elements and running masses 


(of Table |) 


used are obtained from ref. 




VcKM 


/ 0.9747- 0.9759 
0.218 - 0.224 
V 0.003 - 0.018 


0.218-0.224 
0.9735 - 0.9751 
0.030 -0.054 


0.002 - 0.005 \ 
0.032 - 0.054 
0.9985 - 0.9995 ) 


(27) 


Lepton masses of course have no dependence on as to one loop, and their rj values 
are tabulated in Table 


"'^We would like to point out that a more sensible convention to take would be to extract 
lOGeV). This could avoid threshold ambiguities when running up to frit and recent studies 
suggest that renormalons introduce an intrinsic ambiguity into what one means by the pole mass of 
the lighter quarks at these scales anyway. 
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rrit / GeV 

Ve 

Vi. 

Vr 

140 

1.018 

1.018 

1.016 

170 

1.019 

1.019 

1.017 

200 

1.020 

1.020 

1.018 


Table 6 : rji of the Leptons at various rrit. 



alpha3{MZ) 


Figure 3: Renormalisable 3rd family Yukawa coupling for various values of 

rhb 

4.2 The Third Family: Yukawa Unification|5 

The third family have the largest and only renormalisable Yukawa term in this scheme, 
which looks like 

O 33 = (28) 

where A 33 is the universal Yukawa coupling at Mx, unifying 

Xt{Mx) = Xb{Mx) = Xr{Mx) = A 33 . (29) 

Mx is taken to be 10^® GeV. In fact, the results turn out to be insensitive to whether 
we choose Mx = 10^® or 10^^ GeV, for example. 

The third family Yukawa couplings are run from fht to Mx- This is achieved with 

®This subject has been widely considered in the literature (cf. refs. H, H, 1^). We discuss it here 
for completeness. 
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Figure 4: tan/3 prediction for various values of mb. 


the following SUSY one loop renormalisation group (RG) equations (cf. ref. 0): 



(30) 


which are valid in the MSSM, assumed to be the correct theory between fht and Mx- 
We ignore low energy threshold corrections, which should be smaller than the other 
theoretical uncertainties involved (these are briefly presented in section 5.) 


The procedure we follow for the third family is very similar to reh^] for bottom- 


tau Yukawa unihcation, but are here extended to include full top-bottom-tau Yukawa 
unihcation, according to Eq.^. Briefly, the idea is to input frib, fhr and as{Mz) and 
then to predict fht and tan (3 as outputs using the constraint of Yukawa unihcation, 
as in Eq.^ after running the 3rd family Yukawa couplings up to Mx- In practice 
this is complicated by the fact that the Yukawa couplings Xtifht), \b{mt), Xr(fht) all 
depend upon the input value of fht and tan (3. Thus we pick reasonable estimates of 
these quantities to input into the RG routine. The output values obtained from this 
are fed back as inputs until the iteration converges. In this way we select fht, tan/3 
consistent for Yukawa unihcation consistent with given values of fhb, fhr, as{Mz). 
This results in the predictions for mi(pole), Xssi^Mx) and tan/3 shown in Figs. |^-§. 

As Figs. illustrate, the results are highly dependent upon as{Mz), and fairly 
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alpha_s (M_Z) 


Figure 5: If for the third family, as dehned in eq.|3^. The lower (upper) triplet of 
curves is associated with rhb = 4.1 (4.4) GeV. 


dependent on ffib- The value of tan/? is high (35 to about 65) and mt(pole) ranges 
from 130 to 190 GeVQ. Where the curves on the graphs stop for high as{Mz), one 
of the couplings has become too high and so the model is not valid in this region 
of parameter space. Note that for certain superparticle spectrums, these results can 
be perturbed because the determination of fhbijfib) does not include certain one loop 
Yukawa corrections [^. This can have the effect of lowering the prediction of rrit by 
about 30 GeV, which is still compatible with the GDF result for higher as- 


4.3 First and Second family: Diagonal Yukawa Couplings 


In dealing with the hrst and second families we have to confront the problem that the 
Yukawa matrices are not diagonal. As discussed widely elsewhere mE , it is most 
convenient to diagonalise the Yukawa matrices at Mx before running them down to 


rrit- It is then possible to obtain RG equations for both the diagonalYukawa couplings 


An,c,t, ^d,s,b, and the Gabibbo-Kobayashi-Maskawa (GKM) matrix elements | VjlQ 


(ref.0). At one-loop these RG equations can be numerically integrated so that the 
low energy physical couplings have a simple scaling behaviour : 


®This is consistent with the CDF measurement of mt in ref, [pll . 

^The empirical values of \Vij\ were taken to be at fht instead of Mz, introducing an error whose 
magnitude is always less than 1 percent for our analysis. 
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Figure 6: Empirically derived ratios of first and second to third family lepton Yukawa 
couplings at Mx for ffih = 4.25 GeV. 


'A 


U,C 


'A 




^d,s 


mt 


^d,s 


Afi 


Ah 


'A 


e,/! 


'A 


e,/! 


A. 


mt 


A. 


Mx 


Mx 


Mx 




^3Ib+It 


Mr 


I \mx ^ ^1,,+It 
I ^cb I fht 

with identical scaling behaviour to Vet of 14 b, Vts, Vtd, where 


4 = 


pin Mx 


/ In mt 


' Xt (ty 

dvr 


dt, 


(31) 


(32) 


and f = Inp, /i being the MS scale. To a consistent level of approximation 14<i, 14^, 
14s) Ved, Vtb, Xu/Xc, Xd/Xs and Ag/A^ are RG invariant. The GP violating quantity J 
scales as 14|. The I integrals of the third family are displayed in Fig.^for the allowed 
range of ffit and as{Mz). 

Using Eqs. |2^ - p6| we may determine the diagonal Yukawa couplings in this model at 
ffit- The third family Yukawa couplings at Mx (all equal) are given by Fig.^ The hrst 
and second family diagonal Yukawa couplings at Mx are then easily obtained from the 
scaling relations in Eq. 0 . using the Ij integrals in Fig.^. These GUT scale Yukawa 
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Figure 7: Empirical bounds on ratios of diagonal Yukawa couplings of the second 
family at Mx and rhb = 4.25 GeV. 

couplings expressed as ratios are shown in Figs. ||,|^,||. The relative magnitude of the 
diagonal Yukawa couplings between the hrst two families and the third is displayed in 


Fig.p, where the ratio should be 0(e) for the hrst family and O(e^) for the second, if 


the assumption of suppression of the mass scales is to be correct. As seen in the hgure, 
this identihes log(e) ~ 0(—1.5). The second family couplings at Mx are shown in 
Fig.|^ and show that : Ac ~ 6 : 20 : 1. The hrst family couplings are displayed 

in Fig.|| and give A^ : Ac : A^ ~ 3 : 1 : 1/30. These are the patterns that must be 
replicated by our model if it going to be phenomenologically viable. In Fig.^we also 
show the absolute value of \Vcb\ at Mx, calculated by running the value at Mz in 
Eq.|23 up to Mx using EqjM|. 

4.4 The Lower 2 by 2 Block of The Yukawa Matrix: The 
Prediction of the Strange Quark Mass 

The heavy 2 by 2 sectors of the Yukawa matrices are considered hrst, in isolation 
from the rest of the Yukawa matrix. This is possible because we shall assume that 
the Yukawa matrices at Mx are all of the form 



YU,D,E ^ 

V 0 0(e) 
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0 

0 (e) 

0 ( 1 ) 


(33) 
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alpha3(M_Z) 


Figure 8: Empirical bounds on ratios of diagonal Yukawa couplings of the first family 
at Mx using rhb = 4.25 GeV. 

where e << 1 and some of the elements may have approximate or exact texture zeroes 
in them. Assuming the form Eq.^ allows us to consider the lower 2 by 2 block of 
the Yukawa matrices hrst. In diagonalising the lower 2 by 2 block separately, we 
introduce corrections of order and so the procedure is consistent to hrst order in e. 

In searching for the minimum number of operators that ht phenomenological con¬ 
straints (so that maximum predictivity is attained) there are some general arguments 
that yield a lower bound on the number of operators. It is clear that there must be 
at least 3 operators in the lower 2 by 2 block since it should have no zero rows or 
columns, or its determinant would be zero, implying that it contains a zero eigen¬ 
value and therefore a zero mass. Also, for a non-zero \Vcb\-, we require there to be an 
operator in the 32 position in the down and/or up matrices. In fact we shall require 
4 operators in the lower 2 by 2 block since the Clebsch coefficients listed in Table ^ 
are not sufficient to describe successfully all the features of Fig.^ in particular the 
charm-muon splitting (A^/Ac ~ 18)|. With 3 operators, excluding the third family 
discussed previously, we have three Yukawa coupling parameters, not including the 
33 coupling which is already calculated. There are four observables connected with 
these parameters, iWfel, and m^. We use iWfel,ffic,as inputs but the only 

prediction we can make is for the strange mass rhs- 

®The discussion of complex phases is left to later. Here we assume that all of the Yukawa 
parameters are real, which happens to make no difference to our analysis of the lower 2 by 2 block 
(see Section 4.4). 
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Figure 9: \Vcb\ calculated at Mx from the low energy value. 


Given the constraints listed above, the only ansatze that successfully satisfy them 
are listed here: 


^2 

^3 

^4 

^5 

^6 

Aj 

^8 


0 ^ 

- 0 ^ 

0 

0 ^ 

033 

0 

0 ^- 

0 ^ 

0 ^ 

O33 

0 

0^- 

0 ^ 

0® 

033 


0 

0 ^ 

0 ^ 

- 0 ^ 

033 


0 

0 ^ 

0 ^ 

- 0 ^ 

033 

qK 

0 ^ ' 


qM 

033 


qK 

0 ^ ' 



033 



0 

0 ^ 

0 ^ 

- 0 ^ 

033 


(34) 

(35) 

(36) 

(37) 

(38) 

(39) 

(40) 

(41) 


Note that can be used in the 32 position of Ai instead of , yielding exactly 
the same results. Also, ansatze replacing with in A^, with O® in A^, 
with or O'^ in Aq_^ and with in Aq^ also yield the same results. A 7 may 
be transposed to yield the same resnlts, as may A^ if is not in the 23 position 
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Figure 10: fhsil GeV) predictions as a result of Eq.^, to be compared with the 
experimental value quoted in Table 

before transposing, as this would predict zero \Vub\- 

From the above ansatze in Eqs.^^-^, the ratio of muon to strange Yukawa cou¬ 
plings at Mx is found to be: 

='■ (^ 2 ) 

\'^s/ Mx 

where I is a ratio of Clebsch coefficients, predicted to be / = 3 (as in the GJ ansatz) 
or / = 4 (a new prediction). 

Ansatze Ai_6 predict / = 3 in Eq.^ and ansatze ^ 7^8 predict I = 4. 

Eq.^ gives a prediction of the quantity Xs{Mx), which can in turn be run down 
to 1 GeV, yielding a prediction for Fig.|T0| shows the predictions for rhs{l GeV), 
which were found to have negligible dependence upon ffib- Note that the I = 4 curve 
in Fig.|T^ corresponding to ansatze Ay s works just as well as the I = 3 GJ prediction. 

4.5 The Upper 2 by 2 Block of The Yukawa Matrix: The 
Prediction of the Down Quark Mass. 

Having diagonalised the lower 2 by 2 block, we have an effective Yukawa coupling 
in the 22 position for each of the Yukawa matrices. We must now add additional 
operators to account for rrie, m^, |14b| and |Fus|- We shall assume that these 
additional operators are n = 2 ones, as discussed previously together with even smaller 


22 



n = 3 operators which we denote generically Naively the minimum number 

of additional operators in the upper 2 by 2 block is 2: one in the 21 position to 
account for |14s| and one on the top row to avoid a massless first family. However, 
it is impossible to account for {Xu/K)mx 1/30 and {Xu/\d)Mx 1/100 using only 
two operators, since the magnitude of the Clebsch coefficients in Table ^ are simply 
not big enough. To circumvent this problem we shall use the operators in Table ^ 
which give a zero contribution to the up mass and a non zero contribution to the down 
and electron mass, namely In order to achieve a non-zero \Vub\ we require an 

operator in the 21 position which gives a non-zero contribution to the up-matrix. To 
provide a phenomenologically viable \Vub\ it turns out (see later) that the Clebsch 
of the down Yukawa coupling in the 21 position has to be 3 times that of the up 
Yukawa coupling in the 21 position. This implies that the 21 operator must be one of 
QDd Qj, If the 12 operator is simply then this predicts a massless up 

quark. In order to obtain a small up quark mass we must add a small third operator 
(denoted to the 11 or 12 positions®. This leads to the successful upper 2 by 2 

block ansatze shown below. With three extra operators in the upper block we have 
three Yukawa coupling parameters, connected with the 5 observables \ Vus\,mu, rrid, me 
and |I4b|- ’iTT'u, and |I4s| are used as inputs and m^, \ Vub\ are predicted^. 

The possible ansatze in the down sector that account for correct \Vus\-,mu-,m(i,me 
and \Vuh\ and also generate CP violation are: 


B2 

Bs 

B, 

Bb 

Be 


0 

QAd X 

0 0^=^ + 
QAd X 

0 

QAd X 

0 0^=3 + 

QDd X 

0 0^=3 + 02 

QDd X 

0 0^=3 + o^ 

QDd X 


(43) 

(44) 

(45) 

(46) 

(47) 

(48) 


®In all of the ansatze Bi-s, the operator can be placed in the 11 position yielding identical 

results. 

^*^However, there appears an unremovable phase that gives CP violation, so the prediction of \ Vub\ 
depends upon this phase (this is explained more comprehensively in section 4.6). 
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k 

Si 

B2 

Ss 

S 4 

S 5 


St 

Ss 

Ai-Q 

4 

16/3 

2 

3 

4 

(3/2) 

(3/2) 

2 

At ,8 

16/3 

(64/9) 

8/3 

4 

16/3 

2 

2 

8/3 


Table 7: k values predicted by the ansatze Ai to when combined with Bi 
Note that the bracketed entries predict ffidil GeV) to be outside the empirical 


and so are not included in Fig. 11 


to i? 8 - 
range, 



alpha3{MZ) 


Figure 11: md(l GeV) prediction as a result of the Eq.^ with the k values outlined 
in Table labeling each prediction. 


Bs 


0 

QMd X 

0 

QMd X 


(49) 

(50) 


X stands for whatever is left in the 22 position, after the lower 2 by 2 submatrix 
has been diagonalised. Each of the successful ansatze gives a prediction for the down 
Yukawa coupling at Mx in terms of the electron Yukawa coupling: 

y] = 

Mx 

where fc = 3 is the Georgi-Jarlskog prediction of ffid- Other viable possibilities found 
in our analysis are k = 2,4, |, y as shown in Table |^. The k values as dehned in 
Eq.|^ depend upon the I value in the heavy submatrix and are displayed in Table ^ 

Fig.^ shows the hve possible predictions for fhd in the separate ansatze. Again, 
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the results were found to have negligible dependence on rrib because this will only 
change the RG running by affecting rht. The large tan/? dependence factors out since 


we relate the mass of the down quarks to the mass of the leptons only, and these have 


the same dependence upon tan/3. It should be made clear that up until now all of 
the Yukawa couplings are assumed to be real and positive, but in order to do a full 
analysis of the CKM matrix including the CP violating phase we shall now drop this 
assumption. 

4.6 The Full 3 by 3 Yukawa Matrix: The Prediction of \Vub\ 
in Terms of the CP-Violating Phase 

So far we have discussed the lower 2 by 2 block and the upper 2 by 2 block of the 
Yukawa matrices, assuming them to be of the form shown in Eq.^, and taking all 
of the couplings to be real and positive. In fact the effective Yukawa couplings must 
be regarded as complex parameters with relative phases between them. This does 
not affect our analysis of third family Yukawa unification. Nor does it affect the 
predictions of the strange and down quark masses which follow from Eqs.|^,^, since 
I and k are simply ratios of Clebsch coefficients. However the precise list of successful 
lower 2 by 2 block ansatze Ai in section 4.4, and the upper 2 by 2 block ansatze Bi 
in section 4.5 was in fact based on a full analysis of the 3 by 3 Yukawa matrices, 
including complex phases. We shall not repeat the full analysis which led to the 
successful ansatze Ai and Bi here. However in order to illustrate our approach and 
discuss the remaining CKM parameters, it is instructive to consider one example as 
described below. 

The successful ansatze consist of any of the lower 2 by 2 blocks Ai combined with 
any of the upper 2 by 2 blocks Bi, subject to the restrictions shown in Table |^. For 
example let us consider Ai in the lower 2 by 2 block combined with any of the Bi in 
the upper 2 by 2 block, focusing particularly on Ai combined with Bi. Just above 
Mx, before the H, H helds develop VEVs, we have the operators 



(52) 


which implies that at Mx the Yukawa matrices are of the form 



0 7/32X32 ’ e*'?^32 i733e*‘^33 


(53) 
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where we have factored out the phases of the operators and Hik are the magnitudes of 
the coupling constant associated with Oik- Note the real Clebsch coefficients xfk^’^ 
give the splittings between In our particular case Ai, Bi they are given by 


rrU 

■^12 

= 0 

rf.D _ 
"''12 ~ 

1 

- 

"''12 ~ 

"''21 

= 1 

"''21 ~ 

3 

x^ - 
"''21 ~ 

"''22 

= 1 

"''22 ~ 

-1 

x’^ - 
"''22 ~ 

lU 
^ 22 

= 1 

/D 

^ 22 — 

1 

/E 

^ 22 — 

"''32 

= 1 

rfD - 

"''32 ~ 

-1 

J.E _ 

"''32 ~ 


( 54 ) 


and x '^2 7^ 0- make the transformation in the 22 element of Eqj5^ 

xgH22e‘*” - 


( 56 ) 


where H 22 ) 4>22 are real positive parameters. It follows from the Clebsch struc¬ 
ture in Eq.^ that H 22 = ^H §2 4^22 — 4 > 22 - general we shall write H §2 = IH§ 2 , 

where / = 3 in this case. 


At Mx, we have the freedom to rotate the phases of F* and Fj, since this leaves 
the lagrangian of the high energy theory invariant. In doing this we rotate away 5 
phases in the matrices since there are only 5 relative phases: 


■ Fi ■ 


g - h032 -< Al2 ) g 

0 ■ 


F2 


Q Q - K ' t > 32 -< t > 22 ) 

0 


. Fz _ 


1 

0 

0 

1 _ 


■ ■ 


Q - i {- 4 > 32 +< t § 2 ~< t 2 l ) g 

0 

- 

F2 


0 e*'^32 

0 


Fs 


1 

0 

0 

gifliss 



Fi ' 

F 2 
F 3 . 

' Fi 
F 2 
F3 


( 56 ) 


Below Mx, the multiplets F, F are no longer connected by the gauge symmetry in 
the effective held theory, since it is The Standard Model. We now dehne our notation 
as regards the effective held theory below Mx as follows. The ehective quark Yukawa 
terms are written (suppressing all indices) 

{URrY^QLh 2 + {DnYY^QLh, + H.c. (57) 


We transform to the quark mass basis by introducing four 3 by 3 unitary matrices 
^Ul r then the Yukawa terms become 

{URrV^^Vu^Y^V^^Vu,QLh 2 + {DRYV^^Vn^Y^V^^VD.QLh + H.c. (58) 
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where and = Ynj^Y^V})^ are the diagonalised Yukawa 

matrices. With the dehnitions in Eqs.^,Q the CKM matrix is of the form 




CKM 


= 


( 59 ) 


In all of the cases considered, Xi 2 = 0 and x'^’^ = 0 so that the Yukawa matrices 
which result from Eqs.^,55,5(; are 


0 

H 12 XY 2 

0 



H 21 XYI 

Hi 

0 



0 

Hz2xi2 

H 33 _ 



0 

Hi2xf2 

0 



H 21 X 21 

IHi 

0 



0 

H32xi2 

H 33 . 



0 

Hl2x'i2 

gi(0l2'- 

<Pl2) 

0 

X^,H21 

HgY 

{4>22~^2: 

1 ) 

0 


0 H32X^2 -^33 . 


(60) 


In order to diagonalise the quark Yukawa matrices, we hrst make real, by multi¬ 
plying by phase matrices 


g *012 

0 

0 ■ 


g *‘/’22 

0 

0 ■ 

0 

g *</<22 

0 

YU 

0 

1 

0 

0 

0 

1 _ 


0 

0 

1 _ 


(61) 


where we have dehned 022 = 0^2 ~ ^22 ^12 = <P'i 2 ~ 4’i2- This amounts to a phase 

redehnition of the {UjiY and Ul helds. 

To diagonalise the real matrices obtained from the above phase rotations, we hrst 
diagonalise the heavy 2 by 2 submatrices, then the light submatrices as shown below. 


YD 

YU 


0 

(M 

1_ 


0 

0 

-S2 C2 0 


0 64 S4 

_ 0 0 1 


-1 

ICO 

1 

0 

_1 

1 

0 

1_ 


1 

0 

0 

_1 

0 

1 


0 

CO 

Coi 

CO 

0 0 1 


1 

0 

1 

CO 

CO 

1_ 


-1 

0 

0 


C2 -S2 0 

0 C4 —S4 


S2 C2 0 

0 S4 C4 


0 0 1 _ 

- 1 

0 

0 

_1 


Cl —Si 0 

0 C3 -S3 


Si Cl 0 

0 S3 C3 


0 0 1 


( 62 ) 


where q, s* refer to cos6i and sin0j respectively. Note that since sym¬ 

metric Cj, Si are independent of Cj, Sj. 


The diagonal Yukawa couplings of the strange quark and muon obtained from 
Eq.|^ are {Xs)mx = = IH §2 since the 22 eigenvalues are just the 22 
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Figure 12: \Vub\ bounds predicted in terms of the complex phase 0 for various as{Mz) 
and the ansatze ^i-e, corresponding to / = 3. 


elements in this case. The hrst family diagonal Yukawa couplings for the down quark 
and electron are related by 



1 =l 

Mx 


rytD /y>C) 
^21^12 
E E ' 
^21^12 


(63) 


We identify the right hand side of Eq.|^ with k in Eq.|^. The angles are given by S 4 = 

-X32H32/H33, S2 = -X 2 iH 2 i/{Xs)mx, Si = -X 2 iH 2 i/{Xc)mx a^d S3 = -x^ 2 H 32 /H 33 . 

Note that {Xu)mx = —x'X 2 ^ 2 iHi 2 H 2 i/{Xc)mx is small and in the limit O 12 —*■ 0 the 
up quark is massless in the model. 


Denoting 6^ = 6^ — 6*4 and 0 = —022 , we substitute the diagonalising matrices 
from Eqs. ^ and ^ into the CKM matrix in Eq.^ to obtain 


VcKM 


C2Cie*‘^ + S2S1C3 —S2Cie*'^ + S1C2C3 S1S3 
-siC2e*'^ + C1S2C3 siS2e*'^ + C2C3C1 S3C1 
-S2S3 -C2S3 C3 


(64) 


Note that tan^i = This is a generic feature of all the ansatze. To obtain a 
prediction for \Vub\ we note that 




^ \ Vub{mt)\ ~ 


— S2Cie*'^ + S1C2C3 
Vub 11 

Feb 2'21 '^s 



\Vus{mt)\\Vcb{mt)\ 



xO\,(Mx) 

x^^\s{Mx) 


(65) 
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Figure 13: \Vub\ bounds predicted in terms of the complex phase 0 for various as{Mz) 
and the ansatze ^ 7 , 8 , corresponding to / = 4. 


Eq-IHI predicts a value of \Vub\ that is dependent upon the value of 1. This is be¬ 
cause of the appearance of As in Eq.^, which is predicted to have different values in 
Eq.^ depending on 1. The \Vub\ predicted is displayed in Figs. and and hts 
the phenomenological values of 0.002-0.005 successfully only for large values of the 


complex phase. Clearly Figs. p!^ , p!3| predict large values of \Vub\- We emphasise that 
this prediction applies to all of the successful ansatze Ai,Bi. 


5 Conclusions 


We have discussed the problem of fermion masses in the supersymmetric SU(4)® 
SU( 2 )l®SU( 2 ) 7 j model, where the gauge group is assumed to be broken to the stan¬ 
dard gauge group at Mx ~ 10^® GeV. Although the gauge group is not unihed at Mx, 
it is hoped that the model may be embedded in some string theory near the Planck 
scale. Since the model involves no adjoint representations of the gauge group, and 
has no doublet-triplet splitting problem, the prospects for achieving string unihcation 
in this model are very good, and some attempts in this direction have already been 
made []^. However here we have restricted ourselves to the low-energy effective held 
theory near the scale Mx ~ 10^® GeV, and parameterised the effects of string unih¬ 
cation by non-renormalisable operators whose coefficients are suppressed by powers 
of (Mx/M), where M > Mx is some higher scale associated with string physics. 
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We have assumed that the heavy third family receives its mass from a single 
renormalisable Yukawa coupling in the superpotential. The model predicts third 
family Yukawa unihcation (i.e. the top, bottom and tau Yukawa couplings are all 
equal at Mx) leading to predictions for the top mass and tan/3. The other Yukawa 
couplings result from the effect of non-renormalisable operators of order 
where n = 1 operators are suitable for the lower 2 by 2 block of the Yukawa matrices, 
and n = 2, 3 are suitable for the upper 2 by 2 block. In fact we have classihed all 
possible operators in this model for n = 1 and n = 2 . 


The successful ansatz At, Bi in Eqs.|3J-^,^^^ involve 8 real parameters i/ 33 , H 32 , 
H 22 ^, H 21 , Hi 2 , Hi 2 , plus an unremovable phase. With these 8 parameters we can 
describe the 13 physical quantities fhd, ffig, fhc, fhh, m*, rhe, rhr, di, 2 , 3 , 0- 
Third family Yukawa unihcation led to a prediction for mj(pole) = 130 — 200 GeV 
and tan/3 = 35 — 65, depending on as{Mz) and rh},. More accurate predictions could 
be obtained if the error on as{Mz) and rhb were reduced. The analysis of the lower 
2 by 2 block led to 2 possible predictions for rhg depending on whether / = 3 or 4, as 
shown in Fig. 0 (/ = 3 is the GJ prediction). In the upper 2 by 2 block analysis we 
were led to 5 possible predictions for fhd, depending on whether k = 2, 8/3, 3,4,16/3, 
as shown in Fig.|^ (again A; = 3 is the GJ prediction.) Finally, we have 2 predictions 
for \Vuh\ depending on whether Z = 3 or 4, as shown in Figs.[l^Jl^. Both predict a 
large value of | 14 b|, depending on <p. 


The high values of tan/3 required by our model (also predicted in SO(IO)) can 
be arranged by a suitable choice of soft SUSY breaking parameters as discussed in 
ref. , although this leads to a moderate hne tuning problem . The high value of 
tan /3 is not stable under radiative corrections unless some other mechanism such as 
extra approximate symmetries are invoked, may have been overestimated, since 
for high tan /3, the equations for the running of the Yukawa couplings in the MSSM 
can get corrections of a signihcant size from Higgsino-stop and gluino-sbottom loops. 
The size of this effect depends upon the mass spectrum and may be as much as 30 
GeV. For our results to be quantitatively correct, the sparticle corrections to rrif, 
must be small. This could happen in a scenario with non-universal soft parameters, 
for example. Not included in our analysis are threshold effects, at low or high energies. 
These could alter our results by several per cent and so it should be borne in mind 
that all of the mass predictions have a signihcant uncertainty in them. It is also 
unclear how reliable 3 loop perturbative QGD at 1 GeV is. 


Gompared to SO(IO) the lack of predictivity of our model is somewhat dis- 
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couraging. In the SO(IO) model, the spectrum is described by just 4 operators, 
whereas in our model the spectrum is described by 7 operators. One basic reason for 
this is that, unlike SO(IO), our Yukawa matrices are inherently asymmetric. In order 
to £11 out our Yukawa matrices, we need to add operators in the ij and ji positions 
separately. This of course permits asymmetric texture zeroes such as those in Ai_8, 
which have not been studied before. In the upper 2 by 2 block, the SO(IO) model can 
satisfy all of the phenomenological constraints with just 2 operators, the up Yukawa 
coupling becoming very small through a small Clebsch ratio (1/27)^. This permits 
the SO(IO) model to make predictions for the up quark mass and the complex phase 
as well as \Vub\- Whereas we cannot predict in this model, a natural explanation 
is given for its relatively small value in terms of a higher dimensional operator. Note 
that simply not having this operator would not alter any of the predictions, except 
that the up quark would be massless, thus solving the strong CP problem. Thus a 
simple and natural way of obtaining a massless up quark is given by our model with 
Qn= 2 , g gqg which would reduce the number of operators in our model 

by one. 


Despite the lack of predictivity of the model compared to SO(IO), the SU(4) 
®SU(2)£,(8)SU(2 )j:j model has the twin advantages of having no doublet-triplet split¬ 
ting problem, and containing no adjoint representations, making the model technically 
simpler to embed into a realistic string theory. Although both these problems can be 
addressed in the SO (10) model |2^, we find it encouraging that such problems 
do not arise in the first place in the SU(4)®SU(2)i®SU(2)R model. Of course there 
are other models which also share these advantages such as hipped SU(5) or even 
the standard model. However, at the field theory level, such models do not lead to 
Yukawa unification, or have precise Clebsch relations between the operators describ¬ 
ing the light fermion masses. It is the combination of all of the attractive features 
mentioned above which singles out the present model for serious consideration. 
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Appendix 1 : n = 1 Operators 


Following is a list of all n = 1 non-renormalisable operators. These operators are 
constructed from various group theoretic contractions of the following hve helds, 



( 66 ) 


It is useful to dehne some SU(4) invariant tensors C, and SU(2)/j invariant tensors R 
as follows: 



jTQ: 



(67) 


where S^, eap^x^ Sy, e^z are the usual invariant tensors of SU(4), SU(2)/j. The SU(4) 
indices on (Ti, 6 , 10,15 are contracted with the SU(4) indices on two helds to combine 
them into 1, 6, 10, 15 representations of SU(4) respectively. Similarly, the SU(2)/^ 
indices on R 13 are contracted with SU(2)i? indices on two of the helds to combine 
them into 1, 3 representation of SU(2)/j. 

The operators in Tables are then given explicitly by contracting Eq.^ with 
the invariant tensors of Eq.^ in the following manner: 
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0-' ~ (c,)^(Ctn(Ji3)^(Ji3)Z0ZS 

o'” ~ (C,5)S(Ci5)*(fii)i(ii0;O^»” (68) 

Appendix 2: n = 2 Operators 


The n = 2 operators are formed from the following seven helds: 


r)Pip2P3aoq _ rppiaTT jp TTP20 tt TJPiQh^ 
^7l7273mpn ~ -'-'7im-f 72P-'-' -'-'73)1-'-' "-q- 


(69) 


Apart from the invariants dehned in Eq.^ we also define 


(C33,r^z 

( TD \Tnnt 

l-^4jfecfc 


+ syZsj + s'XK + sz€sj + s“y,s2 + s-j^sz 


SiTnsinsit I Sim siTi Sit i emcnet i sm^rn^rt , Sim. s:n Sit \ Simcnsit 

Oc + Oc Ok\ + 0^ dfcA + Ok Oc\ + Ok ObO^ + df, d^dk, 


(70) 


where the possible SU(4) structures which multiply Eq.^ are 


(Avzjli = (Ci5)S2(Ci5)sf(c.)); 
(b)SSm = (Ci5)2?(Ci5)sf(Ci)):; 
(C)’;s;» = (Ci5)22(Ci5)3f(c.)2 

= (Ci)2(Ci)2(Ci)3 

(■E)223 = (Ci)’;(Ci)2:>(c,)g 

1-^^P1P2P3 “ l‘^20jpip2P3l^20ja/37 

('/^')717273 _ ('/^ '\a/3 , (/^ \7172 ,^^73^0- 

v^/pip2p3 v^l0/pjp2^ap3/^o'V^10/ci;/3 ^ 

('tr '1717273 _ '|7172 iJ 73 pcr 

/plp2p3 \^10jp3p2^Qplpcrl^^l0ja;/3 ^ 

/'r'l7l7273 _ '17172 0173 per 

Jplp2p3 l^l0jpjp3^Q:p2pO'l,^10ja;/3 ^ 

z' T'\7l7273 _ '\o/5 , 'i7173^iiZ 72 per 

V'^/piP2P3 l^l0jpjp2*^“P3pe’'V^10/ljl/3 ^ 

|'/y'l717273 _ le /'(e-l^^'l7l73 iJ72po- 

l'^/piP2P3 l^l0jp3p2^opipO'V^10/(jj/3 0 

i'r'\717273 _ '\a/3 , (/^ '17173 , 1 ^ 72 per 

V'^/plp2p3 V^10/pjp3^Q:p2pO'l,^10/ci;/3 ^ 

(3i)2;S3 = (Ci„)“f„(c.„)y(C02 
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f/\/'')7l7273 
/P1P2P3 

= (c'lo: 

vO/9 '>7172 

' P3P2 v^IO/ck/3 

f/n'1717273 
v'^/plp2p3 

= (c^io: 

v«/9 \7372 

hiP2 v^lOJo/3 

/ 717273 

V-' /plp2p3 

= (c^io: 

vO/9 'i7372 

' P3P2 V^10Jq!/3 

f/0'\7i7273 

vVypip2p3 


( p\7l7273 

V /PlP2P3 


( O')717273 
\*^/plp2p3 


/yi\ 717273 
/plp2p3 



mil 

73 

P3 

73 

Pi 

72 

P3 

72 

pi' 


The SU(2)ii; structures which multiply Eq.|^ are 


(«)™7 

= 

f TP \mnt( T 3 \bck^sp^ 
\^i)bck y^AJsqr^ ^to 

(&)r7 

= 

(TP \bck(jp \pns mt 
\^4:)toq\^^)hck^ ^sr 

( \mnp 
y^Joqr 

= 

(jp \bck(Tp \pmn 
\-^‘i-)oqr\^‘^)bck 

idr7 

= 

( jp \mns(jp \bck tp 

y^ijbck 

( \mnp 
y^Joqr 

= 

( rp \bck(jp \pns tm 
{^i)tor\^4)bck^ ^sq 

/ £\mnp 
\J Joqr 

= 



= 


ih)Tq7 

= 

(R.r^e^Pero 

(■\mnp 

y^/oqr 

= 

iRi)7^%o 

( j\mnp 
\J ) oqr 

= 


(U\mnp 
\^/ oqr 

= 

{R^rjiR^rjiRiYr 

n\mnp 

\^)oqr 

= 

iRsfoKRsrJmq 


= 

{Rsyj{R3)TsiRir 


= 

{R,ry{R,)fyR(y 

(oK/ 

= 

iR,rMliRs)7e,,e^- 

ip)o7 

= 

{Ry-{Rs)%{R3)7e,oe^- 

iq)o7 

= 

{RPl{Rsy7R3rs7,oe^^ 

( \mnp 
\ Joqr 

= 

{R^y{Rs)7{R3)Zy-esq 

( \mnp 
Joqr 

= 

{R,y{R,r7{Rs)Zy^esq 

/j.\mnp 
\^)oqr 

= 

{RyiiRsy^RsY^y^e,,. 


The resulting 400 n = 2 operators are of the form 

qES _ (z:)'n'y2'r3 (X\mpnQpip2p3aoq 

^ \^Jpip2p3\^Jaoq '^7i7273mp725 
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where E = A, B .. .T and 6 = a,b.. .t. 
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